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We analyze the scaling behavior at and near a quantum critical point separating a semimetallic from a super- 
fluid phase. To this end we compute the renormalization group flow for a model of attractively interacting 
electrons with a linear dispersion around a single Dirac point. We study both ground state and finite temper- 
ature properties. In two dimensions, the electrons and the order parameter fluctuations exhibit power-law 
scaling with anomalous scaling dimensions. The quasi-particle weight and the Fermi velocity vanish at the 
quantum critical point. The order parameter correlation length turns out to be infinite everywhere in the 
semimetallic ground state. 

Copyright line will be provided by the publisher 



1 Introduction 

Quantum phase transitions in interacting electron systems are traditionally described by an effective order 
parameter theory, which was pioneered by Hertz [|T| and Millis (2). In that approach, an order parameter 
field is introduced via a Hubbard-Stratonovich transformation and the electrons are subsequently integrated 
out. The resulting effective action for the order parameter is then truncated at quartic order and analyzed 
by standard scaling and renormalization group (RG) techniques. 

However, more recent studies revealed that the Hertz-Millis approach is often not applicable, especially 
in low dimensional systems [3]|4]]. For electron systems with a Fermi surface the electronic excitation 
spectrum is gapless. As a consequence, integrating out the electrons may lead to singular interactions in 
the effective order parameter action, which cannot be approximated by a local quartic term. Therefore it 
is better to keep the electronic degrees of freedom in the theory, treating them on equal footing with the 
bosonic order parameter field. Several coupled boson-fermion systems exhibiting quantum criticality have 
been analyzed in the last decade by various methods fS HTOl . 

Recently, a Dirac cone model describing attractively interacting electrons with a linear energy-momen- 
tum dispersion was introduced to model a continuous quantum phase transition from a semimetal to a 
superfluid JTT). The scaling behavior at the quantum critical point (QCP) was studied by coupled boson- 
fermion flow equations derived within the functional RG framework. It was shown that electrons and 
bosons acquire anomalous scaling dimensions in dimensions d < 3, implying non-Fermi liquid behavior 
and non-Gaussian order parameter fluctuations. 

In this work we extend the analysis of the Dirac cone model in various directions, with a focus on 
the two-dimensional case. First, we allow for a renormalization of the Fermi velocity of the electrons, 
which was omitted in Ref. ifTTIl . but indeed turns out to be important. Second, we study the behavior upon 
approaching the QCP from the semimetallic phase at zero temperature. While the pairing susceptibility 
exhibits the expected power-law scaling, we find that the correlation length is infinite everywhere in the 
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semimetallic phase. Finally, we compute the scaling behavior of the susceptibility and the correlation 
length in the finite temperature quantum critical regime. 

The paper is structured as follows. In Sec. 2 we define the Dirac cone model and the corresponding 
action. The derivation of the flow equations is described in Sec. 3, and results are presented in Sec. 4. We 
conclude with a short summary in Sec. 5. 



2 Dirac cone model 

We consider a model of electrons with a linear dispersion relation e^ a = au/|k|, with a = ±1, cor- 
responding to two "Dirac cones" with positive (a = 1) and negative (a = — 1) energy. The chemical 
potential is chosen as /i = 0, such that in the absence of interactions states with negative energy are filled, 
while states with positive energy are empty. The Fermi surface thus consists of only one point, the "Dirac 
point" at k = 0, where the two Dirac cones touch. The action of the interacting system with a local 
attractive interaction U < is given by ifTTI 



S[tp,1p] = I Ipkouri-iko + eiaJlpkaa + U / lp~k,ailpk+ q , a ^k> +q,a>^-k> ,a> I 

r J ka J h' ol 1 J q 



where ip and ^ are fermionic fields. The variables k = (ko, k) and q = (qo, q) collect Matsubara frequen- 
cies and momenta, and we use the short-hand notation j k = T ^ fco j ^J^ d for momentum integrals and 
frequency sums; J k includes also the sum over the band index a and J, includes in addition the spin 
sum over a = f, I . Momentum integrations are restricted by the ultraviolet cutoff Vf |k| < Aq. 

In Ref. ifTTIl it was tacitly assumed that the interaction does not shift the upper and lower Dirac cone 
with respect to each other. To compensate for self-energy contributions which in fact do generate such a 
shift, we have added a fermionic mass term with a U -dependent mass m a to the action S[ip, xjj]. This term 
is tuned such that the Dirac cones touch each other at k = for any U. 

The kinetic energy in Eq. (JT|i is a toy version of the dispersion for electrons moving on a honeycomb 
lattice as in graphene, where the momentum dependence is entangled with a pseudospin degree of freedom 
related to the two-atom structure of the unit cell lfP21 . Note that the kinetic energy and the interaction in 
Eq. (Q} are both diagonal in the spin indices. By contrast, in Dirac fermion models describing surface states 
of certain three-dimensional topological insulators the spin orientation is correlated with the momentum 
lfT3lfT4l . We are not aware of a physical realization of the model Eq. (Q]i in a real material. The model 
was designed to analyze the quantum phase transition between a semimetal and a superfluid in the simplest 
possible setting. Although the model ([T} is reminiscent of the Gross-Neveu model |[T5l . it is not equivalent 
to it. In particular, for the Gross-Neveu model there is no choice of a spinor basis in which the kinetic and 
potential energies are both spin-diagonal. 

The attractive interaction favors spin singlet pairing IfTTI . Therefore, we decouple the interaction in 
the s-wave spin-singlet pairing channel by introducing a complex bosonic Hubbard-Stratonovich field (j> 
conjugate to the bilinear composite of fermionic fields U J, V'fc+</,atV'-fc.c4 ■ This yields a functional 
integral over ip, tp and <fi with the fermion-boson action 

S[4>,^,<j)} = j tpkaai-iko + 6k« + m a ) tpkaa ~ [ 4* q JJ (j) q 
Jkaa Jq U 

($-k,oc\.i>k+qM $1 + ^k+q.a^-k,al 4>* q ) , 

(2) 



ka J q 



where (j)* is the complex conjugate of </>. The boson mass 8 = —\/U > Ois the control parameter for 
the quantum phase transition. In mean-field theory a continuous transition between the semimetallic and 
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a superfluid phase occurs at the quantum critical point t/^! F = — 27rv 2 ./Ao in two dimensions IfTTII . For 
technical reasons explained in Sec. 3, we will supplement the bosonic part of the action by adding a term 
of the form J cp* (Z^qg + Ai,q 2 )(f> q , which regularizes the flow at high scales without influencing the low- 
energy properties of the system. The extra term corresponds to a replacement of the local interaction U 
by a g-dependent interaction U(q) = U/[l — XJiZ^q^ + A^q 2 )), which decreases at large momenta and 
frequencies. From now on we set Vf = 1. 



3 Renormalization group 

Our aim is to derive scaling properties of the electrons and the order parameter fluctuations near the 
quantum phase transition. To this end we derive flow equations for the scale-dependent effective action 
F A [tp, ■tp, </>] within the functional RG framework for fermionic and bosonic degrees of freedom |[T6l4T9l . 
Starting from the bare fermion-boson action r A=A ° [ip, ?/?, <f>\ = S[ip, ip, <p] in Eq. ©, fermionic and bosonic 
fluctuations are integrated simultaneously, proceeding from higher to lower scales as parametrized by the 
continuous flow parameter A. In the infrared limit A — » 0, the fully renormalized effective action T[ip, %p, cp) 
is obtained. The flow of T A is governed by the exact functional flow equation lfl6l 

where r' 2 ) A denotes the second functional derivative with respect to the fields and R A is the infrared 
regulator (to be specified below). The supertrace (Str) traces over all indices, with an additional minus sign 
for fermionic contractions. 



3.1 Truncation 

The functional flow equation Eq. (0 cannot be solved exactly. We therefore truncate the effective action 
with the objective to capture the essential renormalization effects. Our ansatz for T A is a slight generaliza- 
tion of the truncation used in Ref. ifTTI of the following form 



+ r^.^ + rj^| 4 + r A 2( £* , (4) 



where 



1^ = / 4' ka<J (-iZfk + Afe ka +m^)i' kaa , (5) 

J kaa 



= / 0* {Z£q£ + AW + <5 A ) <p q , (6) 

Jq 

U A f 

r |0|4 = ~T <t>q+p4>*q>-p<t>q'4>q , (7) 

" J q,q' ,p 

= g A / / {-Ip-k^k+qM $q + 1pk+q,acfll>-kM <P* q ) ■ (8) 
J hex J q 

The momentum and frequency dependence of r A »^, and also the bosonic interaction rAu, are generated by 
fermionic fluctuations. The fermion-boson vertex r^ 2 0* is actually not renormalized within our truncation. 
The usual one-loop vertex correction, which is formally of order g 3 , vanishes in the normal phase due to 
particle conservation |[T9l . Hence, the coupling g remains invariant at its bare value g = 1 in the course of 
the flow. 

In Ref. IfTTII a restricted version of the ansatz Eq. (|4} with A A = Z A and A A = Z A was used, since it was 
assumed that frequency and momentum dependences renormalize similarly. However, a closer inspection 
reveals that this is not the case. In particular, it turns out that one-loop contributions to the flow of A A 

Copyright line will be provided by the publisher 



4 



B. Obert, S. Takei, and W. Metzner: Anomalous criticality near semimetal-to-superfluid transition 



cancel, while Z A flows to infinity at the QCP. This asymmetry between momentum and frequency scaling 
generates also a significant difference between A A and Z A . 

The initial conditions for the fermionic renormalization factors are Zf° = Af° = 1. The initial 
condition for the bosonic mass is <5 A ° = —1/U, and the quartic bosonic interaction u is initially zero. 
The initial conditions for Z b and A b corresponding to the bare action in Eq. (f2) are 

z b° = A b" = °- 

However, starting the flow with Z^° = A^ a = leads to very large transient anomalous dimensions at the 
initial stage of the flow (for A near An), which complicates the analysis in a (high energy) regime which 
is physically not interesting. The qualitative behavior of the low energy flow (A <C An) and the critical 
exponents do not depend on the initial values of Z\, and A\>. We therefore add a term J </>q(<Z 2 + q 2 ) 4>q 

to the bare action, corresponding to initial values Z^ = A A ° = 1. This term regularizes the model by 
suppressing the interaction for large momentum and energy transfers. 

As regulators in the flow equation ((3) we choose momentum dependent Litim functions l20l . supple- 
mented by a mass shift for the fermions, 

Rj a (k) = A f [-Asgn(e kQ ) + e kct ] 6>(A - |e kQ |) + <5m A , 
itf(q) = A(-A 2 + q 2 )tf(A 2 -q 2 ) , (9) 

where 5m A is chosen such that it cancels to a in Eq. (0 at each scale A. Note that we have set Vf = 1, 
such that A is a common momentum cutoff for fermions and bosons. Adding the regulator functions to 
the quadratic terms in the effective action T A yields the inverse of the regularized propagators, which thus 
have the form 

ajM = [izp„ - 4 eto - < + - 1 = , z}ki _ A , aen( ^ )lmx(AA ^ n ■ do) 

of (.) - [-zni - aw - ^ + ntm-> - - z!;qS+Atm L {A ^ )+s , ■ (») 

Symmetry breaking in interacting Fermi systems is often studied by extending the model to an arbitrary 
number of fermion flavors Nf, and expanding in the parameter 1/Nf. Our truncation captures the leading 
contributions for large Nf. The low energy behavior is captured correctly also to leading order in e, where 
e = 3 — d is the deviation from the critical spatial dimension d c = 3, below which anomalous scaling sets 
in. 



3.2 Flow equations 

The flow equations are obtained by inserting the ansatz Eq. dU) for T A into the exact functional flow 
equation Eq. (O and comparing coefficients. For a concise formulation, we use the following short-hand 
notation for a cutoff derivative and loop integration: 

feo J (2tt) s=bJ 



The scale-derivatives of the regulators read 

$.00 = ~ A f 



d A Rj a (k) = -A A sgn(e kQ )0(A-|e kc 



9 A i? A (q) = -2A A A0(A 2 -q 2 ) , (13) 

where terms proportional to d\ A^ and d\A^ are neglected (as usual, see lfl6l ). The contribution from the 
mass shift, d^dm 1 ^, is also discarded. It is formally of higher order (in a loop expansion) than the terms 
kept, and it does not affect the qualitative behavior. Note that the cutoff derivative in Eq. (fT2l acts only on 
the explicit cutoff dependence introduced via the regulator functions. 
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We thus obtain the following equations for the flow of parameters in our ansatz for T 



A. 



KZf = (g A ) 2 l -^G%(q-k)G$(q) 



(14) 



fc=0 



d A A A = 0, (15) 
d A 6 A = (g A ) 2 



f Gj a (k)Gj a (^k) + ^-fG A ( q ), (16) 

Jka z Jq 



= (g A f f Gj a (k + q)G A a (-k) 

Z °% Jka 

= (9 A f f G A a (k + q)G A a (-k) 



2dqf 



ko 



9=0 



9=0 



(17) 
(18) 

d A u A = -4( 5 A ) 4 / [ G A a {~k)] 2 [G A a {k)] 2 + \{u A f f[G A ( q )} 2 , (19) 

J ko. ^ J q 

d A g A = 0. (20) 

The flow equations for Zf, S, Zf,, u, and g are the same as in Ref. Hill . The momentum derivative in the 
flow equation for is with respect to the first (or any other) component of q. All frequency sums and 
momentum integrations in the above flow equations can be performed analytically, both at zero and finite 
temperature. 

Explicit A-dependences in the flow equations can be absorbed by using rescaled variables 

g A = 9 u A = ■ (21) 



At T > one also has to use rescaled temperatures T fc A = T ^ Zb J Ab and T A = T^f^ to absorb A. 
Anomalous dimensions are defined as usual by logarithmic derivatives of the renormalization factors 

a_ d\ogA£ z _ dlogZt A _ dlogAj dhgZf 
Vb ' rflogA ' Vb ' rflogA ' ]f ~ rflogA ' Vf ~ rflogA ' K ' 

Note that 77^ = 0, since A A does not flow. 

4 Results 

We now discuss the scaling behavior as obtained from a solution of the flow equations, focussing mostly 
on the two-dimensional case. Anomalous scaling dimensions occur in dimensions d < 3 IfTTI . We first 
discuss the ground state, including the quantum critical point, and then finite temperatures. Numerical 
results depending on the ultraviolet cutoff Ao will be presented for the choice Ao = 1. 

4. 1 Quantum critical point 

To reach the quantum critical point one has to tune the bare interacting to a special value U qc such that 
the bosonic mass S A scales to zero for A — > 0. In two dimensions we find U qc w —15.646 for Ao = 1, 
which is about a factor 2.5 larger than the mean-field value. For U = U qc the rescaled variables defined in 
Eq. ( f2TT > scale to a non-Gaussian fixed point, with finite anomalous dimensions in any dimension d < 3. 
Since g A does not flow at all, the scale invariance of g A at the fixed point leads to a simple relation between 
the anomalous dimensions, 

% A + nf + vf = 3 - d . (23) 
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1.5 - 




-Log A 

Fig. 1 Flow of anomalous dimensions as a function of A in the semimetallic ground state for a choice of U close to 
the QCP. 



Furthermore, since the flow of Z\, and Ab is determined entirely by a convolution of two fermionic prop- 
agators, the differences of anomalous dimensions for frequency and momentum scaling of fermions and 
bosons are linked by a simple condition, which can be expressed as 

r?f - ^ = 2(Vf ~ Vf) ■ (24) 

Due to r)f = the above relations reduce to ijf + yj? = 3 — d and — rj^ = 2i]f ■ Solving the fixed 
point equations we obtain the numerical values rj^ pa 0.75, rj^ pa 1.25, and ps 0.25 in two dimensions. 
Hence, at the quantum critical point the order parameter exhibits non-Gaussian critical fluctuations with 
different anomalous scaling dimensions for momentum and frequency dependences. Furthermore, the 
fermionic quasiparticle weight (oc Zj ) vanishes, which implies non-Fermi liquid behavior. Since Af 
remains finite, the Fermi velocity also vanishes at the quantum critical point. This last point was missed in 
Ref. [ 1 Q. Due to the different anomalous dimensions for momentum and frequency scaling, the dynamical 
exponent z acquires an anomalous dimension, too. In the bare action S one has Zq = 1 for bosons and 
fermions. At the quantum critical point, we find 

Z A 

zj = 1 + vf ~ Vf = *b = 1 + Vb ~ % ~ 1-25 . (25) 
The equality between 2& and Zf follows from Eq. (l24l >. 



4.2 Semimetallic ground state 

For \ U\ < \U qc \, the bosonic mass S A saturates at a finite value for A — >• 0, corresponding to a finite pairing 
susceptibility \ — Hitla— yo (5 A ) — 1 . The fermionic Z-factor also saturates, such that — > 0. Hence, 
fermionic quasiparticles survive in the semimetallic state. However, A A and Z A do not saturate, but rather 
diverge as A -1 , such that 77^, r\f — > 1. This is illustrated in Fig. 1, where the anomalous dimensions are 
plotted as a function of A for a choice of U close to the QCP. The QCP scaling is seen at intermediate 
scales, before the anomalous dimensions saturate at the asymptotic values rj^ = and 7/^ = r]^ = 1 for 
A — > 0. A finite anomalous dimension away from the critical point is surprising at first sight. However, it 
can be explained quite easily. An explicit calculation shows that the leading small momentum and small 
frequency dependence of the fermionic particle-particle bubble is linear in two dimensions, as long as 
the propagators have a finite quasi-particle weight. In presence of an infrared cutoff this linear behavior 
is replaced by a quadratic behavior (as in our ansatz), but the prefactors of the quadratic terms diverge 
linearly in the limit A — > 0, reflecting thus the true asymptotic behavior. 
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Fig. 2 Temperature dependence of the pairing susceptibility \ an d the correlation length £ in a double-logarithmic 
plot for U = Uqc- At low temperatures the calculated points lie on straight lines, corresponding to power-law behavior. 

The divergences of A£ and Zfr imply that the correlation length and correlation time of pairing fluctu- 
ations are always infinite in the semimetallic ground state, not only at the QCP. This is consistent with the 
observation that the linear momentum and frequency dependence of the particle-particle bubble leads to a 
power-law decay of its Fourier transform at long space or time distances, instead of the usual exponential 
decay. The divergent correlation length suggests that the entire semimetallic phase is in some sense "quan- 
tum critical". This point of view has indeed been adopted in theories of interaction effects in graphene, 
where the particle-hole symmetric (Dirac) point is interpreted as a QCP separating the electron-doped 
from the hole-doped Fermi liquid. Scaling concepts could then be used to compute thermodynamic ETI 
and transport 11221 properties of graphene near the Dirac point. 

The pairing susceptibility \ is generically finite in the semimetallic ground state and diverges upon 
approaching the QCP. From a numerical solution of the flow equations in two dimensions we have obtained 
the power-law 

X (U) cx (\U qc \ - \U\)-™ , with 7o ps 1.725 . (26) 
4.3 Temperature scaling 

We now present results for U = U qc and T > 0, that is, we approach the QCP as a function of temperature. 
At finite temperature the fermionic propagator is cut off by temperature itself, since fermionic Matsubara 
frequencies are bounded by |fco| > ttT, and the bosonic propagator is regularized by the finite bosonic 
mass 5. Hence, the flow of all unsealed variables saturates for A — > 0. Power-laws are obtained for these 
saturated variables as a function of temperature. In particular, 

A b cx T"'^ A , Z b cx , Z f cx T^f , (27) 

with rfc ps 0.60, ?7jf ~ 1.00, and fjj ps 0.20 in two dimensions. 

In Fig. 2 we show the temperature dependence of the susceptibility \ an d the correlation length £, as 
obtained from a numerical solution of the flow equations at various temperatures in two dimensions. The 
susceptibility is given by the inverse bosonic mass S at the end of the flow (A — > 0), the correlation length 
by £ = v Zb/5. Both quantities obey power-laws at low temperatures, namely 

X(T) cx T"T , with 7 w 1.00 , (28) 
£(T) cx T- v , with v ps 0.80 . (29) 
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Note that we use the letters 7 and v for the exponents by applying the classical definition near a thermal 
phase transition \ oc (T — T c )~ 7 and £ oc (T — T c )~ v to the present situation where T c = 0. The 
correlation length exponent obeys v = z^ 1 , which corresponds to a T^ 1 scaling of the correlation time £ T 
in accordance with general scaling arguments for quantum phase transitions. The exponents 7 and v obey 
the classical scaling relation 7 = (2 — r\^)v. 

5 Conclusion 

We have analyzed the critical properties near a quantum phase transition between a semimetallic and a 
superfluid phase in a two-dimensional model of attractively interacting electrons with a Dirac cone dis- 
persion, correcting and extending a previous work IfPTl . We have studied coupled flow equations for the 
fermionic degrees of freedom and the bosonic fluctuations associated with the superfluid order parameter. 
Both fermions and bosons acquire anomalous scaling dimensions at the QCP, corresponding to non-Fermi 
liquid behavior and non-Gaussian pairing fluctuations. Allowing for distinct renormalization factors for 
momentum and frequency scaling, we have found that they differ substantially at the QCP. In particular, 
the Fermi velocity vanishes. We have also analyzed the semimetallic ground state away from the QCP in 
more detail than previously, finding that the correlation length for pairing fluctuations is always infinite, 
not only at the QCP. Finally, we have studied the scaling behavior upon approaching the QCP as a function 
of temperature. The susceptibility and the correlation length obey power-laws in temperature, as expected, 
and the corresponding critical exponents obey the classical scaling relation. 
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